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We report a series of calculations testing the predictions of the Tran-Blaha functional for the 
electronic structure and magnetic properties of condensed systems. We find a general improvement 
in the properties of semiconducting and insulating systems, relative to calculations with standard 
generalized gradient approximations, although this is not always by the same mechanism as other 
approaches such as the quasiparticle GW method. In ZnO the valence bands are narrowed and the 
band gap is increased to a value in much better agreement with experiment. The Zn d states do not 
move to higher binding energy as they do in LDA+U calculations. The functional is effective for 
systems with hydride anions, where correcting self- interaction errors in the Is state is important. 
Similarly, it correctly opens semiconducting gaps in the alkaline earth hexaborides. It correctly 
stabilizes an antiferromagnetic insulating ground state for the undoped cuprate parent CaCu02, 
but seriously degrades the agreement with experiment for ferromagnetic Gd relative to the standard 
local spin density approximation and generalized gradient approximations. This is due to positioning 
of the minority spin 4/ states at too low an energy. Conversly, the position of the La 4/ conduction 
bands of La203 is in reasonable accord with experiment as it is with standard functionals. The 
functional narrows the Fe d bands of the parent compound LaFeAsO of the iron high temperature 
superconductors, while maintaining the high Fe spectral weight near the Fermi energy. 

PACS numbers: 



I. INTRODUCTION 

Density functional theory (DFT) is formally a the- 
ory for the ground state and not for the spectroscopic 
properties, at least not with the interpretation of Kohn- 
Sham eigenvalues as a one electron spectrum. In fact, 
DFT based methods, using modern generalized gradient 
approximations^ ^ (GGA)s give extremely useful total en- 
ergies and structural properties for a wide range of mate- 
rials, but is known to strongly underestimate band gaps 
in most semiconductors and insulators. It also suffers 
from incorrect predictions of metallic or very small band 
gap semiconducting states in Mott insulators, and in the 
case of undoped cuprates (e.g. La2Cu04 and CaCu02) 
produces an incorrect non-magnetic ground state in con- 
trast to the actual antiferromagnetic insulating state. 

An important aspect of standard DFT calculations is 
that they are computationally efficient. This allows de- 
tailed calculations for wide variety of complex materi- 
als including their full realistic chemical and structural 
details. This greatly facilitates comparison with exper- 
iment. Similarly and in spite of the known problems 
with standard DFT calculations, such calculations have 
proved to be very useful in understanding electronic spec- 
tra of solids, both through interpretation of experiments 
and predictions. Methods that retain the computational 
efficiency of standard DFT, but improve the spectro- 
scopic properties are potentially of great interest. In this 
regard, Tran and Blaha have recently developed a semi- 
local functional^ based on a modification of the Becke- 
Johnson functional.^ Importantly, this is not a hybrid 
functional, and therefore scales like standard DFT rather 
than Hartree-Fock in solids. They showed that this func- 



tional gives very much improved band gaps for a vari- 
ety of insulators, including semiconductors, oxides, rare 
gas solids and lithium halides.^ They also found substan- 
tial gaps in accord with experiment for the 3d transition 
metal monoxides, MnO, FeO and especially the proto- 
typical Mott insulator NiO. Calculations of the optical 
properties for a variety of heavier halides also showed 
much improved agreement with available experimental 
data.^ 

Here we present further tests of this functional empha- 
sizing systems that are difficult to describe for a variety of 
different reasons. We also discuss in more detail the elec- 
tronic structure of ZnO, which was one of the systems 
studied by Tran and Blaha. We find general improve- 
ments in the band gaps over the standard PBE func- 
tional in all materials. We also find an antiferromagnetic 
insulating ground state for the Mott insulator CaCu02. 
However, the electronic structure differs from that pre- 
dicted by other methods and its consistency with exper- 
imental data as regards the distribution of the d spectral 
weight is not clear. The electronic structure of ferromag- 
netic Gd is degraded with respect to the standard PBE 
functional. 



II. APPROACH 

The calculations were performed with the stan- 
dard general potential linearized augmented planewave 
(LAPW) method with local orbitals,^ ^ as implemented 
in the WIEN2k package."^ We use the Perdew Burke 
Ernzerhof (PBE) GGA as a reference, with the excep- 
tion of CaCu02, where we use PBE+U calculations, 
with Coulomb parameter U-J=7 eV and self-interaction 
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correctioEF (SIC) and about mean fielcF^ (AMF) double 
counting corrections, since standard GGAs do not pro- 
duce an antiferromagnetic ground state in that caseP^ 
Except as noted, we used experimental lattice parame- 
ters and relaxed all non-symmetry fixed internal coor- 
dinates with the PBE GGA. Also, except as noted, we 
performed calculations in a scalar relativistic approxima- 
tion for the valence states, while core states were treated 
relativistically. The basis sets and zone samplings were 
well converged, with additional local orbit als used to re- 
lax linearization errors. We used the standard LAPW 
matching as opposed to the faster augmented planewave 
plus local orbital (APW+lo) approach. 



III. SEMICONDUCTORS AND 
NON-MAGNETIC INSULATORS 

We begin with ZnO, which is a well studied material 
with a band gap of 3.44 eV. Standard functionals give 
an unusually large band gap error in this material. Cor- 
recting this error is of considerable importance in under- 
standing the defect physics and semiconducting proper- 
ties of ZnO.^^ Our PBE calculations, which we performed 
with the experimental lattice parameters of a=3.2489 A, 
c=5.2049 A, and internal parameter ii=0.381 yield a band 
gap of 0.80 eV. The band structure is shown in Fig. [T] 
The unusually large underestimate is generally thought 
to be associated with a position of the Zn 3d bands at 
too l ow bi nding energy with respect to the valence band 
g^^gg [14115] ^Yiis case, the resulting overly large level re- 
pulsion between the 3d states and valence bands pushes 
the valence band maximum up, reducing the band gap. 
This view is supported by photoemission measurements, 
which plac e the 3d states ~ 7 - 8 eV below the valence 
band edge,^^^^ although it should be noted that this 
could possibly reflect properties of the surface and not 
the bulk. Considering the importance of this material it 
would be desirable to further check the position of the 
Zn d bands, for example by higher energy photoemission 
or other spectroscopies. 

The ban d ga p error is largely corrected by GW 
calculations.lSnSi These calculations place the Zn d at 
high binding energy ~ 6.5 eV, as compared to ~ -5 eV 
in standard PBE calculations (see Fig. |2|. The error in 
the band gap is also greatly improved in LDA+U cal- 
culations, where the Coulomb parameter is used to shift 
the Zn d orbitals to higher binding energy.'^ Hybrid func- 
tional calculations also greatly improve the band gap and 
shift the d states to higher binding energy.!^ We obtained 
a band gap of 2.65 eV, which is much improved relative 
to the PBE functional, but is still lower than experiment. 
This is almost exactly the same value as that reported 
by Tran and Blaha.'S Our band structure (Figs. [I] and 
[2| shows Zn d bands that are narrowed relative to the 
PBE result, but are not shifted to higher binding energy. 
This may well be the origin of the remaining discrepancy 
between the TB-mBJ and experimental band gaps. 




FIG. 1: Comparison of the PBE and TB-mBJ band struc- 
tures for ZnO. 
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FIG. 2: (color online) Comparison of the PBE (top) and TB- 
mBJ (bottom) electronic densities of states for ZnO. The Zn 
d character is by projection onto the Zn LAPW sphere, radius 
2.05 Bohr. 



We next discuss La203. This high-k dielectric ox- 
ide occurs in a hexagonal (P3ml number 164) struc- 
ture with lattice parameters a=3.943 A, and c=6. 141 A. 
Our calculated internal coordinates are ZLa=0.2468 and 
2:01=0.6458, in good agreement with the experimental 
values of 0.2467 and 0.6470, respectively.^^ La203 has a 
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FIG. 3: Comparison of the PBE and TB-mBJ band struc- 
tures for La2 03. 



FIG. 5: Spin orbit split band structure of ferromagnetic Gd 
with the PBE (left) and TB-mBJ (right) functionals. 
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FIG. 4: (color online) Comparison of the PBE (top) and TB- 
mBJ (bottom) electronic densities of states for La2 03. The 
La / character is by projection onto the La LAPW sphere, 
radius 2.5 Bohr. 



low temperature optical band gap of 5.34 eV-^^The cal- 
culated PBE and TB-mBJ band structures and densities 
of states (DOS) are shown in Figs. |3]and[4| respectively. 
The valence bands are derived from O 2p states. In the 
PBE calculation the conduction bands come from La sd 
states that hybridize with La 4/ states that occur just 



above the band edge, while in the TB-mBJ calculation 
the conduction bands come from the La 4/ states. Ob- 
taining a correct band gap in this compound depends on 
a correct placement of the La 4/ resonance. The calcu- 
lated TB-mBJ band gap is 4.74 eV, which is higher than 
the PBE value of 3.85 eV. It is in much better agree- 
ment with experiment but still ~ 0.6 eV too low. Most 
likely this reflects the position of the La 4/ resonance. 
Interestingly, the position of the La 4/ resonance with 
respect to the O 2p derived valence band edge is almost 
the same in the PBE and TB-mBJ calculations, although 
the 4/ bands are narrower with the latter functional be- 
cause the other conduction bands are shifted to higher 
energy. This behavior is reminiscent of what we found 
for La halides.'^ 



IV. FERROMAGNETIC GD 

We next discuss an open / shell material. Gd is a fer- 
romagnetic metal, with a half filled 4/ shell. This leads 
to a high spin moment, nominally 7 /ub from the 4/ elec- 
trons. The actual experimental moment per atom is 7.63 
/i^, which includes a contribution from polarization of 
the other valence states. ^^ ^^ The value of the moment is 
not sensitive to the treatment of the / electrons provided 
that they remain fully polarize d, as they do in standard 
DFT and LDA+U calculations.^^ On the other hand, 
standard DFT calculations produce a Fermi energy den- 
sity of states, N{Ep)^ that is higher than that inferred 
from specific heat measurements. While there can be 
an enhancement of the specific heat in experiments for a 
variety of reasons, such as elect ron-phonon interactions, 
spin fluctuations, etc.^ a reduction is difficult to under- 
stand. The high N{Ef) in calculations comes about be- 
cause the minority spin / states, while above the Fermi 
energy, are close to it, and hybridize with the valence 
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FIG. 6: (color online) Electronic DOS of ferromagnetic Gd. 
Majority and minority spin projections are shown above and 
below the horizontal axis, respectively. The orbital projec- 
tions are onto the Gd LAPW spheres of radius 3.1 Bohr. 



bands. LDA+U calculations move the minority bands to 
higher energy, reducing the hybridization and improving 
agreement with experiment. In effect, standard DFT cal- 
culations underestimate the exchange splitting of the 4/ 
shell. 

We performed calculations for hep ferromagnetic 
Gd, including spin orbit with the experimental lattice 
parameters, a=3.629 A, c=5.796 A. We aligned the 
magnetization with the c-axis. 

The band structures and DOS are shown in Figs. [5]and 
[6j respectively. As may be seen, the TB-mBJ functional 
keeps the average / position approximately the same as in 
the PBE calculations, but reduces the exchange splitting 
opposite to LDA+U calculations. This brings additional 
/ character to the Fermi energy, reduces the moment, 
and increases N{Ef). On a per atom basis we obtain 
N{Ef)=1.73 eV-^ with the PBE functional and 3.64 
eV"-*^ with the TB-mBJ functional. The inferred experi- 
mental specific heat value is 1.57 eV~^.^^ Our calculated 
PBE spin moment is 7.62 /i^ in agreement with experi- 
ment and prior studies, while the TB-mBJ value is 6.65 
fiB- Therefore, it may be concluded that the TB-mBJ 
functional seriously degrades agreement with experiment 



HYDRIDES 



Hydrides present an interesting challenge to DFT cal- 
culations, particularly anionic hydrides. Standard DFT 
calculations have large errors in their descriptions of 
the one electron H atom because of an incomplete can- 
cellation of the Hartree and exchange correlation self- 
interactions. The trihydride, YH3, which is a trans- 
parent insulator, is incorrectly described as a metal in 
the LDA. This has been attributed to self-interaction on 
the H, the ordinary problems in describing band gaps 
within the LDA, other correlation effects, and complica- 
tions with the crystal structure.'^^"^^" Several calculations 
have shown that a semiconducting gap can be opened in 
GW quasiparticle calculations and also the with the non- 
local weighted density approximation (WD A) to DFT.^ 
The WDA is self-interaction free for one electron sys- 
tems. Unfortunately, while the insulating behavior of 
YH3 is likely an electronic effect and not a consequence 
of structural details,'^ the fact that the exact structure 
is uncertain complicates interpretation of calculated re- 
sults. 

Here we present results for three well studied hydrides, 
specifically LiH, MgH2 and AIH3. LiH occurs in the fee 
NaCl structure. Our calculations were done with the ex- 
perimental lattice parameter for LiD at 83K, a=4.0447 
A.^^ MgH2 is tetragonal, with a rutile structure. We 
used the experimental lattice parameters, a=4.5198 A, 
c=3.025 A.^ The internal H coordinate, u was kept at 
its experimental value of 7x=0. 30478 since the force on H 
was found to be small with this value. AIH3 occurs in 
a hexagonal, R3c structure.!^ We used the experimen- 
tal lattice parameters, a=4.449 A, c=l 1.804 A, and re- 
laxed internal coordinate. The structure has Al on site 
6b (0,0,0) and H on site 18 e (x,0,l/4), with x=0.6319, 
in good accord with the reported experimental value of 
x=0.628 (Ref. gg)). 

Our calculated PBE and TB-mBJ band structures are 
given in Figs. [7|[8]and|9l for LiH, MgHs and AIH3, re- 
spectively. The corresponding DOS are in Figs. [lO] and 
[11] for LiH and AIH3. Our PBE band structures are simi- 
lar to those obtained in prior DFT calculat ions .I^^IEII The 
band gaps are of charge transfer character between pre- 
dominantly H s derived valence bands and predominantly 
metal derived conduction bands. 

We start with LiH. Both functionals predict a direct 
gap at the X point. The calculated band gaps are 3.01 
eV and 5.08 eV for the PBE and TB-mBJ functionals, 
respectively. The experimental band gap of LiH is 5.0 
eV as quoted in Ref. 44, Thus the TB-mBJ is in almost 
perfect agreement with experiment for LiH. MgH2 has an 
indirect gap. The PBE and TB-mBJ values are 3.70 eV 
and 5.74 eV, respectively. The experimental band gap of 
MgH2 is 5.6 eV.^^ Therefore, as in LiH, almost perfect 
agreement with experiment is obtained with the TB-mBJ 
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FIG. 7: Band structure of LiH with the PBE (left) and TB- 
mBJ (right) functionals. 



FIG. 9: Band structure of AIH3 with the PBE (left) and 
TB-mBJ (right) functionals. 





FIG. 8: Band structure of MgHs with the PBE (left) and 
TB-mBJ (right) functionals. 



FIG. 10: (color online) Calculated electronic DOS of LiH with 
the PBE and TB-mBJ functionals. 



functional. We are not aware of an experimental band 
gap for AIH3, but recent GW calculations yield a values 
of 3.54 eV or 4.31 eV depending on the treatment. Our 
values are 2.27 eV and 4.31 eV for the PBE and TB-mBJ 
functionals. In the cases studied the valence bands are 
somewhat narrower with the TB-mBJ functional relative 
to the PBE results. The results for hydrides with the 
TB-mBJ functional are very encouraging. 



VI. CABe AND SRBe 

The hexaborides CaBe and SrBe are modest band gap 
semiconductors .E^l'^ Like the anionic hydrides they have 
charge transfer gaps between valence bands derived from 
a compact orbital where self-interactions may be impor- 
tant and conduction bands with metal atom character. 




FIG. 11: (color online) Calculated electronic DOS of AIH3 
with the PBE and TB-mBJ functionals. 
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In this case the valence bands are derived from B 2p 

states. Also, similar to YH3, standard DFT calculation 

predict a semimetallic state due to a band overlap, while 

both WDA and GW calculations yield small insulating 
gaps.ElHSl 

The nature of the bands is, however, different from the 
hydrides. In the hydrides the valence bands are nom- 
inally from the filled H Is state, while the conduction 
bands are nominally from the metal orbitals, although 
there may be some hybridization. In the hexaborides, 
the bonding of the B network plays a central role. The 
valence bands are derived from bonding combinations of 
B p states, while the conduction bands come from mix- 
tures of non-bonding and antibonding B p states with 
Ca/Sr d states. The center of the metal d states is at ~4 
eV above the band gap, but the hybridization of B p and 
Ca/Sr d orbitals persists down to the conduction band 
minimum. The band gap underestimation could there- 
fore come from incomplete cancellation of self-interaction 
errors in the B p derived valence bands, or alternatively 
an overestimate of the Bp- Ca/Sr d hybridization that 
then pushes the conduction band minimum down (analo- 
gous to the pushing up of the valence band edge in ZnO), 
or a combination of these two possible errors. 

We did calculations for CaBe and SrBe using the ex- 
perimental simple cubic crystal structures. We did not 
relax the internal B coordinate because the forces with 
the experimental values were small. The lattice param- 
eters used were 4.1514 A, for CaBe, and 4.1953 A, for 
SrBe. The internal coordinates were x=0.2019 for CaB^ 
and x=0.2031 for SrBe. Our calculated band structures 
for CaBe and SrBg are given in Figs. [T2]andJT3) respec- 



tively. The corresponding DOS are in Figs. [M^and 15 



We find a metallic state in PBE calculations due to 
band overlap at the X point in both materials. On the 
other hand, with the TB-mBJ functional we obtain small 
gaps: 0.10 eV for CaBe and 0.18 eV for SrBe. There is 
a wide range in measured band gap values for the hexa- 
borides. However, recent optical measurements for CaBe 
yield a band gap of 0.25 eV,^^ while prior measurements 
indicate a larger gap ~1 eV.^^ Considering the T de- 
pendent resistivity, and the spectroscopic measurements, 
the TB-mBJ result is clearly better than that of the PBE 
calculation, but the gap is underestimated. In contrast 
to the hydrides, we do not find a valence band narrowing 
with the TB-mBJ functional. Also, the position of the 
Ca/Sr d states in the DOS (not shown) at ~4 eV is almost 
exactly the same in the PBE and TB-mBJ calculations. 
It will be desirable to check this against spectroscopic 
data, if such data becomes available. An upward shift of 
the d bands would increase the band gap. 



VII. CACUO2 




FIG. 12: Band structure of CaBe 
TB-mBJ (right) functionals. 



with the PBE (left) and 




FIG. 13: Band structure of SrBe 
TB-mBJ (right) functionals. 



with the PBE (left) and 




As mentioned, Tran and Blaha reported calculations 
for the transition metal monoxides, MnO, FeO and NiO 
finding improved band gaps along with changes in the d 



FIG. 14: (color online) Calculated electronic DOS of CaBe 
with the PBE and TB-mBJ functionals. 



7 




FIG. 15: (color online) Calculated electronic DOS of SrBe 
with the PBE and TB-mBJ functionals. 



density of states of NiO that differ from those found in 
LDA+U calculations P Here we present calculations for 
CaCu02. This material is the simplest parent compound 
representative of undoped cuprate superconductors. 

These materials are Mott insulators when undoped, 
and become high temperature superconductors when 
doped. Importantly, they have a non-trivial band struc- 
ture associated with strong hybridization between Cu d 
orbitals, particularly the a bonding Cu dr^i_yi -Op com- 
binations. These combinations lead to an antibonding 
band, of dx2_y2 orbital symmetry crossing the Fermi en- 
ergy and giving rise to the Fermi surface on which super- 
conductivity occurs. This hybridization also leads to a 
large interatomic exchange coupling, and high magnetic 
energy scale, which may be important for the high super- 
conducting critical temperatures. When undoped, this 
band is half filled, and as mentioned the Fermi surface is 
destroyed in favor of a Mott insulator. While the Fermi 
surfaces of the doped phases are generally well described 
in DFT c alculations, the magnetism in these materials is 
not.fni^In contrast to MnO, FeO and NiO, these mate- 
rials are described as non-magnetic in DFT calculations, 
and so the ground state is qualitatively wrong, not only 
in terms of electronic structure but also in the fact that 
there is no moment formation. This can be corrected 
by LDA+U (or GGA+U) methods, which yield an in- 
sulating antiferromagnetic ground state. These meth- 
ods also shift the d spectral weight away from the Fermi 
energy to form what may be termed upper and lower 
Hubbard bands. This shift of d spectral weight away 
from the Fermi energy is consistent with photoemission 
measurements.'^ However, while there is a shift towards 
Hubbard bands at high binding energy in both doped and 
undoped cuprates, there does remain d spectral weight at 
the band edge in the insulators and Fermi energy in the 
superconductors. 

We performed calculations using the PBE+U method, 
with U-J=7 eV and two double counting corrections, the 
self-interaction correction (SIC) scheme and the about 



mean field (AMF) scheme. We used the experimen- 
tal tetragonal crystal structure (there are no internal 
coordinates).!^ For the magnetic order we assumed a 
G-type ground state, i.e. the in-plane nearest neigh- 
bor checkerboard antiferromagnetism that is common to 
cuprates with antiferromagnetic stacking along the c- 
axis. 

The resulting electronic DOS are compared with re- 
sults obtained using the TB-mBJ functional (and no 
added U) in Fig. 
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First of all we note that unlike 
standard DFT treatments the TB-mBJ functional suc- 
cessfully stabilizes the antiferromagnetic insulating state. 
The calculated spin moments as measured by the polar- 
ization inside the Cu LAPW sphere (radius 2.0 Bohr) are 
0.636 /iB, 0.695 //b, and 0.646 j^b for the PBE+U (SIC), 
PBE+U (AMF) and TB-mBJ functionals, respectively. 
The reduction from the ideal spin 1/2 value of 1 /i^ re- 
flects hybridization with the O p states. The calculated 
band gaps are sensitive to the value of U-J assumed in 
the PBE+U calculations. With the value of 7 eV used 
here, we obtained gaps of 1.5 eV and 1.3 eV with the SIC 
and AMF double counting schemes, respectively. The 
TB-mBJ gap is comparable at 1.7 eV. The essential dif- 
ference between the PBE+U and TB-mBJ results is seen 
in the distribution of the Cu d DOS. The two double 
counting schemes, yield somewhat different distributions 
of the d DOS, but both have a strong shift of the oc- 
cupied d DOS to higher binding energy. The TB-mBJ 
DOS instead places substantial d DOS near the valence 
band edge. This is similar to what Tran and Blaha found 
in NiO. In cuprates, it is clear that there are Hubbard 
bands at high binding energy, which are present in the 
PBE+U but not the TB-mBJ calculations. However, it 
is unclear from existing experimental data whether the 
description near the valence band edge is better with the 
PBE+U or the TB-mBJ functionals. One possible av- 
enue for resolving this is to note that when projected 
onto a given site, the Cu d character near the band edge 
with the TB-mBJ functional has a substantial spin po- 
larization, which could in principle be detected in spin 
polarized scanning probe experiments. 



VIII. LAFEASO 

The iron based superconductors provide an interest- 
ing contrast to the cuprates.'^ Unlike cuprates, these 
materials show strong d spectral weight around the 
Fermi energy, in accord with DFT calculations, weak to 
moderate correlations, and no Hubbard bands. DFT 
calculations,^^ show an Fe d^ metal, with a high density 
of states and a semimetallic character, i.e. small Fermi 
surfaces: hole-like sheets around the zone center, and two 
electron sheets at the zone corner. These features are in 
accord with experiment, and are degraded by methods 
that build in strong correlations. One difference between 
DFT calculations and experiment is that the d bands are 
narrower in experiment. 
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FIG. 16: (color online) Electronic DOS and projections for an- 
tiferromagnetic CaCu02 with the PBE+U method (see text) 
and the TB-mBJ functional. The Cu d projection is onto the 
Cu LAPW sphere of radius 2.0 Bohr. Majority and minority 
spin projections are shown above and below the horizontal 
axis, respectively. 



We performed calculations for LaFeAsO using the ex- 
perimental 175 K crystal structure of de la Cruz and 
CO- workers,^ including the experimental internal coordi- 
nates. We did not relax the internal coordinates because 
they are known to be incorrectly given by non-magnetic 
GGA calculations, while magnetic calculations incorpo- 
rating the spin density wave give overly strong mag- 
netism, presumab ly du e to spin fluctuations neglected in 
DFT calculations .l^^l^ In any case, this should be kept 
in mind when comparing calculations, since the band 
structure, especially near the zone center, is sensitive 
to the As position in the unit cell. In particular, there 
are additional hole sheets of Fermi surface that depend 
on whether one uses the experimental or calculated As 




FIG. 18: Band structure of LaFeAsO with the PBE (left) 
and TB-mBJ (right) functionals. 
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FIG. 19: (color online) Fermi surface of LaFeAsO with the 
PBE (left) and TB-mBJ (right) functionals. The shading is 
by band velocity and the corners are F points. 



position. Here we focus on the large scale features of 
the electronic structure that are not dependent on this 
choice. These are the energy dependent character of the 
bands and the shape of the main Fermi surface, which 
governs superconductivity.^^ 

The calculated non-spin-polarized DOS, band struc- 
ture and Fermi surfaces are shown in Figs. T7| 18 and 
[191 respectively. As may be seen, the bands within 2 
eV of Ef are of mainly Fe d character with modest hy- 
bridization of As p states in both the PBE and TB-mBJ 
calculations. The As and O bands are at higher binding 
energy, yielding nominally Fe^+ metallic ^ sheets. Also 
in both calculations there is a pseudogap in the DOS 
starting near corresponding to the semimetallic char- 
acter. This is seen in the small disconnected Fermi sur- 
faces. One difference is that the TB-mBJ calculation has 
a ~10% smaller d band width, which would marginally 
improve agreement with experiment. The detailed dis- 
persions of the bands near Ep are different between the 
two calculations. The electron pockets are larger in the 
TB-mBJ calculation. Also, there is one fewer small hole 
Fermi surface section. However, while these features may 
be important to understanding the materials, they are 



at the level of differences that arise from different treat- 
ments of the As position, and we do emphasize them here. 
The overall agreement of the PBE and TB-mBJ results 
in this material is gratifying. 



IX. SUMMARY AND CONCLUSIONS 

We performed calculations for a diverse set of materials 
using the recently developed semilocal density functional 
of Tran and Blaha. The present results taken in com- 
bination with the tests that they reported^ and calcu- 
lations for halides^ strongly support the conclusion that 
this functional greatly improves the band gaps and elec- 
tronic structure of simple semiconductors and insulators. 
It also greatly improves upon standard density function- 
als in that it stabilizes an antiferromagnetic insulating 
ground state in CaCu02. However, the quality of the 
description of d bands still needs clarification. In ZnO, 
the Zn Zd states may be too shallow, and in CaCu02 
the distribution the d DOS in the valence bands differs 
from that obtained with the standard LDA+U method. 
The large scale features of the TB-mBJ electronic struc- 
ture of LaFeAsO are similar to those with the PBE func- 
tional, although there are changes in details that affect 
the hole Fermi surface sections. Turning to / band ma- 
terials, the properties of La2 03 are improved relative to 
PBE calculations, although the 4/ bands are hardly af- 
fected at all by the change in functional. On the other 
hand, the properties of ferromagnetic Gd are seriously 
degraded with respect to the PBE functional. 

In summary, the results indicate that the TB-mBJ 
functional is a very useful development for treating the 
electronic structure of simple semiconductors and insula- 
tors, including hydrides, while care is needed in applying 
it to d and / band materials. 
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